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Abstract 

Duality principle for approximation of geometrical objects (also 
known as Eudoxus exhaustion method) was extended and perfected 
by Archimedes in his famous tractate "Measurement of circle". The 
main idea of the approximation method by Archimedes is to construct 
a sequence of pairs of inscribed and circumscribed polygons (polyhe- 
dra) which approximate curvilinear convex body. This sequence al- 
lows to approximate length of curve, as well as area and volume of 
the bodies and to obtain error estimates for approximation. In this 
work it is shown that a sequence of pairs of locally polar polyhedra 
allows to construct piecewise-affine approximation to scherical Gauss 
map, to construct convergent pointwise approximations to mean and 
Gauss curvature, as well as to obtain natural discretizations of bending 
energies. 

Keywords: polar polyhedra, discrete curvatures, surface of bounded curvature, 
bending energy. 

Discrete curvature functionals and surfaces of bounded curva- 
ture. One of the hard problems of modern geometry is approximation 
of nonregular surfaces by polyhedra. In the sense of intrinsic metric (based 
on distance along surface) this problem was solved in the works of A.D. 
Alexandrov and his scientific school [Ij. A.D. Alexandrov developed theory 
of "good" approximation of manifolds of bounded curvature by polyhedral 
manifolds. However these results are not sufficient to establish "good" con- 
vergence in the sense of extrinsic metric. The class of surfaces being man- 
ifolds of bounded curvatures in the intrinsic sense is well defined: they are 
called surfaces of bounded curvature[2], [3], [4], [5j. 

Extrinsic curvatures for polyhedra are introduced using integral rela- 
tions. Gauss-Bonnet theorem allows to assign to the vertex of polyhedron 
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curvature which is equal to angular excess of its conical neighborhood [12j, 
[1]. Balance equations for vector mean curvature can be used to derive 
discrete mean curvatures for polyhedra and to construct discrete approxi- 
mation to Laplace-Belt rami operator. To this end one can also use variation 
of surface area and its relation with the sweep volume [9j. In [7] with each 
region on the surface it is associated a tensor which in the smooth case is 
the average of curvature tensor over this region. For polyhedral domain 
the same value provides weakly convergent estimator of the curvature ten- 
sor. In [8] it is shown that if a sequence of polyhedral surfaces converges 
to a regular surface in Hausdorff distance, then the following conditions are 
equivalent: a) convergence of normal fields, b) convergence of metric tensors, 
c) convergence of area, d) convergence of Laplace-Beltrami operators. 

In order to define "good" approximation by polyhedral surface, one have 
to define such curvature measures for non-regular surface M, which can 
be introduced via sequence of polyhedral surfaces Pk which the following 
properties: a) Pk converge to M pointwise when k +cxo; b) Pk converge to 
M uniformly in intrinsic metric; b) positive and negative parts of curvature 
of Pk converge to to positive and negative part of curvature of M in a weak 
sense; c) spherical Gauss map of Pk converge to spherical map of M in a 
weak sense. Approximation problem in such a setting is still not solved [lOj. 

Consider curvature functional which can be used for investigation of 
non-regular surfaces. Let M denote regular 2D surface in (regular in 
a sense that it admits thrice continuously differentiable local parameteriza- 
tion). Consider functional 



where da is the surface area element, A ^ M?^'^ is the matrix of the shape 
operator or curvature tensor defined by equality 



Here G and T are matrices of the first and second fundamental forms of 
the surface, respectively, and g{A) is certain curvature density measure. 
Functional Eg{M) is called bending energy of the surface. If energy Eg{M) 
is bounded, then absolute Gauss curvature of the surface is bounded as well. 

Matrix A is nothing else but the jacobian matrix of the spherical map 
11. Let us remind that spherical map identifies with each point p of regular 
surface M a point h — v{p) on a unit sphere where v{p) is the unit 
normal to surface. One can compute intersection point q between plane, 
passing through h and orthogonal to v{p) and a ray going through 
where point p' belongs to some neighborhood of which is shown in fig. [1] 
a). Mapping p' ^ q \s called normal map and defines normal image of 
neighborhood of p. 
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Fig. [T] a) Spherical and normal images, b) normal graph over the surface. 

It is said that polyhedral surface is normal graph over M if projection 
of surface M onto along normals to M is homemorphism, which is 
shown on fig. [T]b). 

Well known example of bending energy is given by mean total quadratic 
curvature measure 

E2(M) = tT{A^A) da 

M 

Absolute minimum of this functional is attained when a surface homeomor- 
phic to sphere is precisely the sphere. Mean quadratic curvature measure is 
not suitable for description of nonregular surfaces since it is not defined for 
polyhedra. In other words it takes infinite value for polyhedral surface. 

If bending energy majorates absolute curvature and remains bounded for 
refined sequence of polyhedra then one can expect that the limiting surface 
for this sequence will be surface of bounded curvature. 

One can consider the following curvature measures which make sense for 
polyhedra: 

Ei{M) = J{{tT{A^A))^ + \detA\)da (2) 

M 

and 

E,SM) = j.ftriA-A)^i±lj^^. (3) 

v2j^ {e + tT{A^A))2 

where ^ > is a constant. 

Duality principle and approximation of surfaces by polyhedra. Let 

us consider the method for construction of discrete bending energies suitable 
for approximation of nonregular surfaces. 
Consider 2D paraboloid 

P = {x : xs = u{xi,X2),u{xi,X2) = 
= ^{huxj + 2/ii2XiX2 + h22xl)} 
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We shall use upper index / to denote vectors from R^, while values without / 
will denote their orthogonal projections onto plane X3 = 0. It is convenient 
to write function u as u{p) — ^p^Hp^ where H is the shape operator matrix 
of paraboloid P at the origin. 



a) b) 




Fig. [2] a) Polyhedral surface inscribed into elliptic paraboloid; b) dual face and normal 

image of the vertex. 

Consider fragment of convex polyhedral surface inscribed into elliptic 
paraboloid. The faces of this polyhedron incident to vertex p{ = are shown 
in fig. [2] a). The plane X3 = is tangent plane at the origin. Tangent planes 
passing through the vertices p^j lying at the edges, incident to p\, cut on the 
plane X3 = polygon Qi. In fact Qi is the face of dual (polar) polyhedra 
P^ circumscribed around the same paraboloid. Let us denote by V{p\) the 
set of vertices of P^ belonging to edges, incident to p\^ while notation V(G) 
is used for the set of vertices of the face G. 

One can construct normal image of the vertex p\^ namely the convex 
polygon Fi on the plane X3 = 1. The vertices of this polygon are intersections 
of rays passing through p\ and orthogonal to faces, incident to p\^ with the 
plane X3 = 1. Polygons Qi and Fi are shown in fig. [2]b). 

Consider vertex g[ of the polygon Qi. Vector g[ can be found as the 
solution of the linear system 

n\p\f{q{-p\) = 

n\p)f{qi-p]) = 0, j G V(G,), i^j 

where v) are unsealed normals to paraboloid and p^- are vertices of a face 
Gk G star(p^), i.e. incident to p\. 

Linear system @ is nonsingular. In order to obtain normal n\ one can 
compute gradient of the function X3 — u{xi^X2) — 0, i.e., 

n\p]) - ( ) , n'(p^) = ( ; ) 
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Thus linear system (j4]) can be reduced to 

i-Hpjfiqk -Pj) - u{pj) = 0, 

or 

{Hqkfpj = j e V(Gfc), i ^ j (5) 

It can be easily verified that if the number of equations in system ([5j) exceeds 
the number of unknowns it remains consistent. 

Now let us find vertex of the polygon Fi. Its coordinates are the 
solution of linear system 

(p'j - P'd'^ifi - P^) = 0, je V{Gk), i + j 
This system simplifies to 

pjfk + u{pj) = 0, i e V{Gk), i + j (7) 

Prom equations ((S)), ((T]) follows equality 

fk = -Hqk, (8) 

which hold for all vertices of polygon Fi. Thus the following theorem is 
proven: 

Theorem 1. (JT4^I) polygons Fi and Qi are affine equivalent, i.e. Qi = 
(j)l{Fi) and jacobian matrix of affine map coincides with —H, where H 
is the matrix of shape operator of paraboloid P at the origin. 

Let us note that deriving equality fk — —Hq^ we did not used the fact 
that matrix H is positive definite. Formally ([8]) holds for arbitrary matrix 
H . It is just required that matrix with vectors pj as columns has the full 
rank. Thus duality principle for computation of curvature tensor can be 
applied in the case of hyperbolic paraboloid, shown in fig. [3l 

a) b) 




Fig. O a) Polyhedral surface inscribed into hyperbolic paraboloid; b) dual face and 

normal image of the vertex. 

From the duality principle it follows that face Gk of polyhedron P/^, 
incident to p\ corresponds to vertices and If faces Gm and Gk have 
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common edge, then vertices and should be connected by an edge. The 
same is true for vertices and These arguments can be applied in the 
case when boundary of Fi and Qi is self-intersecting closed polyline. In this 
case vertex p[ can be called nonregular since polyhedron provides poor 
local approximation to paraboloid P in the neighborhood of p\. 

The condition that P^ is inscribed polyhedral surface and P^ is the cir- 
cumscribed one is the particular case of polarity with respect to paraboloid 
It is well known that relation of polarity, i.e. one-to-one correspon- 
dence between point and a plane can be introduced using arbitrary surface 
S of the second order [12j. 

Consider rays originating from certain point p and tangent to S. All 
points of contact with S will lie in the same plane which is polar to point 
p. It may happen that point p is situated in such a way that tangent rays 
cannot originate from it. Then one can draw an arbitrary plane 11 through p. 
Tangent rays passing through intersection points between S and 11 consitute 
a cone. When plane 11 is varied, the summit of this cone sweep the plane 
which is precisely the polar plane to point p. Let us remark that in convex 
analysis and optimization only special case of polarity with respect to sphere 
is considered. 

Suppose now that vertices p^j of P^ do not lie on the surface of paraboloid 
P, i.e. 

{p^j)s = ^pjHpj + Sj 
The plane of the face, polar to p^j is defined by equality 

(^')3 + iPj)3 = pjHx 

Suppose that pi = 0, then {p\)s = S^. Point fl is defined by system ([6]), i.e. 

pjfk + IpjHpj + Sj -Si^O, je V{Gk), i + j (9) 

while Qk is defined as the intersection point of a planes, polar to the vertices 
of edge Gk, thus 

{Hq.fpj = ^pjHpj + Sj - Si, j e V(Gfc), i ^ j (10) 

As a result one obtains equality fk — —Hq^^ which means that theorem [T] 
holds for polyhedral surfaces, polar with respect to paraboloid P, and not 
just for inscribed and circumscribed polyhedra. 

Now consider the case when there exists face Gk of P/^, parallel to the 
plane X3 = 0. Vertex q^, dual to this face is intersection of planes polar 
to points p^j being vertices of G/e, and q^ = 0. Consider normal image of 
the neighborhood of vertex g^, i.e. polygon Bj^ with vertices computed as 
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solutions to linear system 

n\GkY{h]-qi) = W{Gk)l 
- Qinfo^ - = 0, ieV(Q^), i^i, ^''^ 

where n\Gk) denotes normal to face Gk- Elementary calculation shows that 
i.e. it is proven the following theorem 

Theorem 2. (1141) Polygons Bj^ and Gk are affine equivalent, i.e. — 
^k{Gk) and jacobian matrix of affine mapping (j)k coincides with —H, where 
H is the curvature tensor of paraboloid P at the origin. 

Let us consider arbitrary regular closed 2D surface M and polyhedral 
surface P^ inscribed into M. One can choose cartesian frame Xi such that 
surface M can be locally written as X3 = f{xi^ X3), and 

/(xi, X2) = u{xi,X2) + 0(|xp), hij = ^^^^ (0, 0), 

which means that P is touching paraboloid at the point p\ of the surface M. 
Denote by Fi normal image of neighborhood of the vertex p\. The vertices of 
polygon Fi are computed via solution of linear system ([6]). Let us construct 
dual polyhedral surface P^, consisting of the faces tangent to M. Tangent 
face, dual to p\ is denoted by Qi. Vertices of polygon Qi can be found by 
solving ([4]), where are normals to the surface M. In local coordinate 
frame they can be written as 

n\p]) = ( -^[^'^^ ) 

For a vertex of a dual polyhedral surface P^ one can find in turn the 
normal image, polygon Bj^. The plane of this polygon is parallel to the face 
G/e, dual to Vertices of polygon B^ can be found using equation ffTTD . 

Definition 1. The vertex p\ of polyhedral surface P^ is called regular if 
its dual polygon Qi and orthogonal projection of polygon Fi onto the plane 
of polygon Qi are star-shaped domains with respect to the point p\. 

Definition 2. The vertex ql. of polyhedral surface P^ is called regular if 
face Gk of Ph and orthogonal projection of polygon Bi onto the plane of face 
Gk are star-shaped domains with respect to the projection of point q^,. 
It is clear that these two definitions are mutually symmetric. 
For non-convex surface the above definitions are too restrictive. One can 
use definition of weak regularity: 
Definition 3. The vertex p\ of polyhedral surface P^ is called weakly 
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regular if its dual polygon Qi and orthogonal projection of polygon Fi onto 
the plane of polygon Qi are simple polygons and contain point p\ inside. If 
polygon Bj^ is simple then one can assume that vertex g[ is weakly regular. 

a) b) 




Fig. m Triangula! ions of elliptic and hyperbolic paraboloids and their projections onto 

horizontal plane. 




Fig. [5] Dual polyhedral surfaces and their projections. 

On fig. m there are shown fragments of triangulated polyhedral surfaces 
Ph inscribed into elliptic and hyperbolic paraboloid and their projection on 
the plane X3 = 0. Dual polyhedral surfaces and their projections are 
shown in fig. [5l In the case of elliptic paraboloid all dual faces are convex 
polygons, while in the hyperbolic case dual faces Qi are quadrilaterals with 
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concave sides (edges), i.e. the turn of edges from the side of Qi is nonpositive. 

It should be noted that dual polyhedral surface in the case of elliptic 
paraboloid is precisely the Voronoi generatrice and projection of its faces 
onto horizontal plane make partitioning of a plane into convex polygons 
being affine image of Voronoi partitioning. 

In general case one cannot state that polygons Qi and Fi are affine 
equivalent. The polygons Gk and Bj^ are not affine equivalent as well. Thus 
one have to construct piecewise affine homemorphisms (j)^ \ Qi ^ Fi and 
<i>k Gk ^ Bj^. Without loss of generality one can assume that 0^ and (^k 
are 2D mappings. If all vertices of Ph and P^ are regular, then polygons 
Qi and Gk can be triangulated simply by connecting their vertices with 
the point with respect to which they are star-shaped. In the case of weak 
regularity one have to construct more general triangulations. Let us denote 
by triangulation of Qi, and by Tj^ - triangulation of Gk- The images of 
these triangulations under maps 0^ and (f)k are triangulations 7^^ and Tj^ . 
The number of triangles can be sharply reduced if all faces Gk are triangles. 
Then mapping (j)i \ Gk ^ Bk can be taken as affine one. 

Denote by —A^^ jacobian matrix of affine map of triangle G 
onto m-th triangle of 7^^, i.e. 

A^m = -V0*|r4, (12) 

Denote by —Akm jacobian matrix of affine map of triangle Tkm ^ onto 
m-th triangle of 7^^, i.e. 

Akm = -^(t^kkkm (13) 

It is obvious that if all vertices of P^ and P^ are regular, and diameters 
of faces of P^ tend to zero, matrices A^^, Akm converge to exact value of the 
curvature tensor A, i.e. 

where '0'*' is the homeomorphism which maps M onto P^ along normals to 
M. 

Let us remark that matrices A'^^ and Akm general are not symmetric. 
Hence in order to compute principal curvatures one have to use singular 
values of these matrices instead of eigenvalues, and use SVD (singular value 
decomposition) of A'^^ and Akm order to compute principal directions. 

In practice exact surface normals are not known. Many methods are 
available for computation of normals at the vertices of polyhedra. One pos- 
sible way to find normals approximately is to require that discrete curvature 
of Ph is as close as possible to discrete curvature of dual surface P^. In or- 
der to make this comparison possible, one has to construct piecewise-affine 
homeomorphism t/j^ P^ ^ Ph- 
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a) b) 




Fig. [6l Construction of piece wise- affine homeomorphism iph Ph ^ Ph- 

On fig. [6] a) it is shown fragment of polyhedral surface, with direction of 
view being orthogonal to the plane of polygon Qi. Projection of the face Gk 
onto the plane of dual face Qi has nonempty intersection with Qi, namely 
polygon D^, shown in fig. [6] b). If all vertices of and are regular 
then domain consists of two triangles from . Preimage of is the 
quadrilateral D belonging to the face Gk also consisting of two triangles 
from Tj^. Thus mapping j/j^ ' ^ D is affine for each triangle from these 
pairs. 

In the case of weak regularity one can construct more general piecewise- 
linear homeomorphism'^/^. On fig. [71 it is shown fragment of saddle surface. 
In this case dual face Qi is nonconvex and it is too restrictive to require that 
Qi is star-shaped with respect to p\. 




Fig. [T] Construction of piecewise- affine homeomorphism iph ■ Ph ^ Ph, non-convex case. 

In this more general case dual faces can be triangulated in such a way 
that all angles are bounded from below and construction of homeomorphism 
should be based on projection and intersection of general triangulations. 

Denote by 0/^ and 0^ piecewise affine mappings coinciding with phik on 
Gk and with (j)i on Qi, respectively. We will use notations V(/)/^ and 
for piecewise constant functions which coincide with jacobian matrices of 
mappings 0/^ and 0^ where the jacobians are well defined. 

As a result the curvature deviation measure can be defined as follows 

sip) = \\vMMp))-'^<Plip)\\, 

where p is a point lying on P^, and function 6{p) is piecewise constant. Here 
1 1 • 1 1 mean Frobenius matrix norm. 

Theorem 3. Consider closed polyhedral surface P^ which is inscribed into 
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regular closed surface M and is normal graph over M realized via homeor- 
phism '0. Suppose that all faces of are triangles with minimal angle 
hounded from below, while length of edges satisfies quasi-uniformity condi- 
tion 

Ch^lj ^h, < C < 1 

Let closed polyhedral surface be normal graph over M realized via home- 
omorphism , and vertices of P^ lie on the faces of P^. Suppose that all 
vertices of Ph and P^ are regular and angles of triangulations and 
are bounded from below, and the following closeness condition for discrete 
curvatures holds 

sup 5 < 0{h) 

i.e. for all k and for every triangle ^ the following inequality holds 

\\AijVi;h-AlJ\^0{h), 

where = iph{Tij), Tij G . Then with h ^ mappings ip and 
converge to identity, and 

\\A{p'^-Al\\^0{h), 
\\A{r-\qi))-Ak^\\i^O{h) 

One can consider average deviation of discrete measures, namely inte- 
grals 

52{Pn.Pl) = / \\V<l^h[^{p))-V<l^im^da (14) 
Ph 

and 

5i(P,,P,^) = |(||V(/),(^b))-V(/)^b)||+ (15) 
Ph 

^deiVPh{^{p))-deiV4i{p)\)dcj, 

where da denote surface differential on P/^. 

If conditions of theorem [3] hold, but instead of maximum norm deviation 
one imposes closeness in a weak sense, i.e. 

52(P/z,P^)^0when/i^0, 

then one can conjecture that the following weak convergence estimate holds 

j \\A{'4,-\p)) + V<t>h{p)fda^Q, (16) 
Ph 

j \\A{r-\p)) + Vct>l{p)\Ua^O 
PI 



11 



If one assumes 



then one can expect the fohowing convergence estimate 




(17) 



Ph 

+ \det A{xp~'^{p)) -det V(f)h{p)\) da 



and 



J(\\A(r-\p))-mip)\\ + 

+ \det A{i;*-^{p)) -det V<j)l{p)\) da 



(18) 



Thus one can expect that the sequence of polyhedral surfaces where dis- 
cretized functional ([2|) is bounded remains in the class of surfaces of bounded 
curvature. 

Nice property of resulting discrete energy is that it provides exact values 
of curvatures for polyhedron, inscribed into sphere. Consider polyhedra 
inscribed and circumscribed around unit d-dimensional sphere in R^+^. On 
fig. [8] it is shown simplest case of a unit circle (d = 1). 



One can easily see that face Qi, dual to vertex pi of inscribed polyhedron 
Ph is congruent to the normal image of vertex pi. Hence shape operator 
matrix is equal to minus identity matrix. It is also obvious that congruence 
property holds in d-dimensional case as well. 

In 2004 A.I. Bobenko [13] introduced discrete Willmore energy W{Ph) 
(called also conformal energy) being invariant to 3D Mebius transforms. 
This discrete energy is supposed to approximate exact Willmore energy 




F. 



> 

VP* 



Fig. [8l For sphere dual face and normal image are congruent. 
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where A is gradient of spherical map, i.e. shape operator matrix. 

In [1^ it was proven, that W{Ph) — 0, when is convex polyhedron 
inscribed into sphere. 

The same property holds for duality-based energy 

E{P*h) = 5](^tr(A^^) -det^)|Q,area(Qi) (19) 

i 

since matrix A on each dual face Qi is equal to — /. 

Main differences with conformal energy are that the discrete duality- 
based energy properties are the same for d-dimensional sphere in R^+^, and 
terms in discrete Willmore energy (IT9]) converge to the same terms in exact 
Willmore energy without any assumptions on polyhedra and P^ beside 
those guaranteeing area convergence. This is not the case for conformal 
energy which converges only on special polyhedra. 

Duality principle allows to discern polyhedral approximations from limit- 
ing polyhedral surfaces. Let us explain this statement using simple example. 

a) b) 



Fig. [9l a) Dual regular polyhedra, b) coinciding polyhedra. 



Figini a) shows cube and octahedron, dual(polar) to each other with 
respect to sphere. Thus one can assume that limiting surface M is sphere or 
another regular surface. Suppose now that Ph is obtained by adding vertices 
at the centroids of cube faces which is shown in fig. [9] b). Cube faces are 
triangulated including centroids. Formal construction of dual surface P^ 
leads to conclusion that it coincides with P^ thus making limiting surface 
M undistinguishable from P^ and P^. One can compute exact value of 
integral ([2|). 

E,{Ph - p;:) - J2 ie|2tg(Y) + E ^^^<Pii m 

e i 
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where first sum is the sum of discrete mean curvatures concentrated at edges, 
LPe denote angle between normals to faces adjacent to edge e, while second 
sum is taken over vertices. 

It is clear that one can construct discrete curvature functionals using 
spherical image of a neighborhood of a vertex instead of normal image. As a 
result functional ([2j) can be written in the way independent on auxiliary val- 
ues, such as normals at the vertices of polyhedron M, which define polygons 
Fi 

E,iPh = P^) = Yl Ie|^<= + E l^^l' (21) 

e i 

where Kj - absolute extrinsic (Gauss) curvature at z-th vertex of polyhedron. 

Thus using spherical image has some advantages but computation of dis- 
crete curvatures requires complicated formulas of spherical trigonometry , 
while discrete curvature functional based on normal image need only com- 
putation of minors of small matrices. Note that functionals fl20l) always 
majorates (l2T]) . This relation is true for poledral surfaces approximating 
regular ones as well. 

Let us remark that class of surfaces with bounded curvature measures 
in the sense of normal image is the subset of lipschitz continuous surfaces. 
It should be noted as well that method for computing discrete curvature 
described above can be applied in the case of surfaces with boundary. More- 
over, the same duality principle can be applied for d-dimensional surfaces in 
R^+^ and allows to introduce extrinsic discrete curvatures in this case. In 
particular, the proofs of basic theorems [T] and [2] does not use the fact that 
rf= 2. 
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a) b) 




Fig. 1101 a) Dual polyhedron with errors, b) projection of surface triangulation and dual 

surface. 

Dicsrete curvature measures for non-regular polyhedra. Poly- 
hedron inscribed into regular surface can contain non-regular vertices even 
in the case when face normals converge to exact normals of the surface with 
refinement of polyhedral surface. An example of non-regular triangulation 
Ph inscribed into elliptic paraboloid 

xs = 2x1 + 

is shown in fig. [TOl One can see from the figure [10] a) that certain "faces" 
of the dual surface are selfintersecting, while triangulation P^ consists of 
well shaped triangles. 

Projections of surface triangulation P^ and dual surface P^ on the plane 
X3 = are shown on fig. [10] b). Here direction of maximal curvature is 
horizontal. It is clear that non-convex edges of P^ lead to dual edges with 
wrong orientation and hence to self-intersection. Thus eliminating "non- 
convex" edges and creation of another edges in resulting quadrilaterals will 
make all vertices of P^ regular, even though may lead to triangles with small 
angles. 
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Fig. \TT\ Schwarz lantern: triangulation of the cylindrical surface. 

Fig. [TT] shows the so-cahed Schwarz lantern, i.e. polyhedron inscribed 
into circular cylinder. This lantern can be constructed by subdividing the 
side surface of the cylinder by planes orthogonal to its axis into m equal 
parts. Into circle which is obtained in each cross section one inscribes regular 
n-sided polygon. This polygon is rotated by vr/n when passing to the next 
cross-section. As a result the side surface of discretized cylinder consists 
of isosceles triangles. When m, n ^ oc this triangulation converge to the 
surface of cylinder pointwisely, but the limit of the sum of the triangle 
areas is 27tR{H^ + tTr^i?^)^, where K is the limit of the ratio n/m? when 
n — > oc, provided that it exists [15]. Here R^H denote radius and height of 
the cylinder. 

The neighborhood of nonregular vertex p\ of a surface triangulation can 
be described as a "fan" , i.e. as a cone with wrinkles. Normal image 11+ of 
the surface of this fan is non-simplyconnected domain with self-intersecting 
booundary, which is shown in fig. [12] a). If cone belongs to a certain 
half-space, one can construct its convex envelope - cone . Normal image 
of convex cone is convex polygon drawn in fig. 1121 a) by bold lines. 
We will call S+ by principal component of the normal image. Principal 
component can be constructed for saddle point as well, which is shown in 

fig. nab). 

In order to construct principal component one need to eliminate edges 
from K~ until one obtains canonical saddle K~ with a normal image being 
quadrilateral with the sides of nonpositive turn. 
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Generic cone with undefined normal image is shown in fig. [T2l c). 

Extraction of principal component essentially means that discrete cur- 
vature is constructed for some other polyhedral surface with the same set of 
vertices. Thus, even for approximation of regular surface, part of the edges 
only prevent "good" approximation. 

It is natural to call principal component of discrete curvature measure 
the discrete curvature measure based upon principal component of the nor- 
mal image. It should be noted that principal component of the discrete 
curvature measure for Schwarz lantern converges to curvature measure of 
exact cylinder even though area converges to wrong value. 




Fig. 1121 Filtering of edges and principal component of normal image. 

When normal image Fi is self-intersecting, its boundary can be decom- 
posed into simple closed arcs bounding simple polygons with different ori- 
entations. Summing up the areas of these polygons with the signs corre- 
sponding to orientations of their boundary contours, one can obtain a value 
which converges to 27r — ^ Oij with diameter of Fi tending to zero. Here Oij 
denote angles incident to the i-th vertex, and 27r — ^ 9ij is intrinsic (Gauss) 
curvature at the i-th vertex. One can compute the sum of absolutes values 
of areas of simple polygons thus obtaining the value 

areaF/ + 5(F0, > 0, 

where Fi^ denotes principal component of the normal image. In order to 
take into account the original shape of polyhedral surface, one can augment 
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principal component of the integral curvature measure, satisfying inequality 
([1]) by a term 

i 

where the sum is taken over all nonregular vertices. 

It may happen that after certain edges are eliminated by filtering proce- 
dure, the cones whose summits are adjacent vertices of polyhedral surface, 
become inconsistent. In this case one assumes that certain faces of the sur- 
face Ph coincide with faces of dual surface P^. In this case normal images 
and dual faces for these singular vertices should be decomposed into sub- 
domains. Some sub-domains should be eliminated while discrete curvatures 
defined on other subdomains allow to approximate one-sided limits of exact 
curvatures. 

Simply put, one can formulate the optimality criteria for polyhedral ap- 
proximations as follows: optimal polyhedral surface should minimize differ- 
ence between intrinsic absolute curvature and extrinsic absolute curvature. 
Since extrinsic absolute curvature cannot be less than intrinsic one, then to 
some extent this optimality means minimization of absolute extrinsic cur- 
vature but unlike |16| . this minimization is applied only at the non-regular 
vertices of polyhedron. When polyhedral surface is regular in the sense of 
definitions [1] and [2] its intrinsic and extrinsic curvatures coincide. 

One can conclude that duality principle allows the possibility to inves- 
tigate convergence of polyhedral approximations of non-regular surfaces. It 
can be applied in the general case of d-dimensional surface in R^+^. 
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